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Abstract. Homogenization theory is used to calculate the macroscopic dielectric constant from the quantum 
microscopic dielectric function in a periodic medium. The method can be used to calculate any macroscopic 
constitutive relation, but it is illustrated here for the case of electrodynamics of matter. The so-called cell 
problem of homogenization theory is solved and an explicit expression is given for the macroscopic dielectric 
constant in a form akin to the Clausius-Mossotti or Lorentz-Lorenz relation. The validity of this expression 
is checked by showing that the macroscopic dielectric constant is causal and has the expected symmetry 
properties, and that the average of the microscopic energy density is the macroscopic one. Finally, the 
general expression is applied to Bloch eigenstates. Finally, the corresponding many-body problem is briefly 
discussed. 
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1 Introduction 

When light falls onto a crystal, the quantum interaction of 
light with matter is represented locally by a microscopic 
dielectric function e(r,r'). To calculate the macroscopic 
(homogeneous) dielectric constant, one considers that the 
charge distribution created by the light polarizes the crys- 
tal which, in turn, reacts by inducing an electric field that 
modifies the charge distribution. 

In quantum chemistry, this reaction field represents 
the infiuence of the solvent on the solute [^|^,^. Other 
descriptions use the related concept of local fields 
which has been described, for disordered media, by a clus- 
ter expansion 

For periodic media, the local field effect was evaluated 
in the early sixties by Adler and Wiser How- 
ever, many textbooks in solid-state physics still identify 
the macroscopic dielectric constant e with (e), the aver- 
age of the microscopic dielectric function e(r,i;0 over a 
unit cell. Even the most cautious authors do 
not go beyond the relation e — eo(eo + 2(e))/(4eo — (e)), 
that goes by the name of Clausius-Mossotti or Lorentz- 
Lorenz (see Refs.|l^,|l^,0 for a history of this relation). 
Since the local field effect can be quite large [Q, its ne- 
glect can probably be attributed to the numerical burden 
of the standard local field formula . 

In this paper, homogenization theory will be used to 
provide various alternative formulas for the calculation of 
the local field effects. Homogenization theory is tailored to 
calculate the macroscopic e from the microscopic e(r, r'). 



The major trick of the method is to expand all fields as 
a series in ascending powers of the ratio of the lattice pa- 
rameter over the wavelength of the external electromag- 
netic field. Physicists have sometimes used such an expan- 
sion JTqI , but mathematicians exploited it extensively and 
turned it into a rigorous tool. Homogenization theory is 
a branch of applied mathematics that started its expan- 
sion in the late seventies to understand the macro- 
scopic properties of composite, porous, disordered, bub- 
bly, fibrous or layered materials. It is now a fully fledged 
theory that has been applied successfully in many 
areas, such as mechanics, acoustics, electrostatics, fluid 
dynamics, statistical physics, numerical analysis, materi- 
als sciences, electromagnetism |2^, petroleum geophysics 
[p4[, shape- memory alloys ||2^ and pile foundation analy- 
sis M. 



This paper starts with an introduction to homogeniza- 
tion theory, then the microscopic Maxwell equations are 
homogenized to yield the macroscopic Maxwell equations 
and the constitutive relations. Since homogenization the- 
ory is not a usual tool of solid state physics, the calcula- 
tions in these sections will be given in detail. The rest of 
the article will follow the elliptic style of normal research 
papers. Various formulas will be given, corresponding to 
different physical situations. Then, several desirable prop- 
erties of the dielectric constant will be derived. It will be 
shown that the macroscopic dielectric constant is causal, 
has the proper symmetry properties and gives the ex- 
pected energy density and Poynting vector. The formalism 
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Fig. 1. Limit of sin(a;/a) + 6 as a ^ 

is adapted to the case of Bloch eigenstates and a band- 
structure formula for e is given. Finally, the much simpler 
case of the many-body dielectric function is briefly dis- 
cussed. 



2 The point of view of homogenization theory 

Ever since the nineteenth century, physicists have homog- 
enized microscopic systems by performing averages over 
distances very small compared to the macroscopic wave- 
length and very large compared to the atomic dimensions. 
Mathematical homogenization started when homogeniza- 
tion was not considered any longer as an averaging op- 
eration, but as a limit process. As an example, consider 
that a microscopic quantity can be represented by the pe- 
riodic function f{x) — sin a; -\- b. To express the fact that 
the oscillation is very fast, mathematicians tried to give a 
meaning to the limit of f{x/a) as a — > 0. In other words, 
what is the limit of a periodic function when its period 
becomes infinitely small ? 

The contact with the physical point of view comes from 
the fact that the limit is the average of f{x) over a pe- 
riod : f{x/a) — > 6 as a — > 0. However, figure 1 shows that 
this limit is somewhat unusual. Physicists generally do not 
bother very much with the various mathematical limits, 
but in the present case, it is necessary to realize that we 
do not deal with a strong limit, but with a weak limit. 
Strong convergence of f{x) to its limit 6, which means 

that, for a fixed L, dx\f{x/a) — 6| ^ as a ^ 0, is 
clearly not realized for oscillating functions (for our ex- 
ample, /q dx\f{x/a) — 6| — > 2L/7r). Oscillating functions 
enjoy only weak convergence, which is defined by the fact 
that, for any smooth function g{x)^ dxg{x)f{x/a) ~* 

b /J" dxg{x) as a ^ 0. In physical terms, the fast oscilla- 
tions of f{x/a) are smoothed away by measuring f{x/a) 
with a device having a finite resolution g{x). 

Strong convergence enjoys many nice properties, for 
instance the product of two strongly convergent functions 



converges to the product of the limits, but no such thing is 
available for weak convergence. In our example, the weak 
limit of p{x/a) is 6^-1-1/2 and not b^ . Since it is clear 
that the average of the product of two functions is not 
the product of the averages of the functions, the reader 
may wonder why it is useful to consider homogenization 
as a limit instead of an average. The advantage of homog- 
enization as a limit is twofold. On the first hand, it leads 
naturally to an asymptotic expansion of the functions in 
terms of a (i.e. f{x/a) ^ b + af^^\x/a) -|- • • so that 
corrections to the average become available. On the other 
hand, it enables us to treat not only functions but also 
differential equations. 

The way mathematical homogenization deals with dif- 
ferential equations can be illustrated by a simple exam- 
ple |2^]. Consider a wire of length L, uniformly charged 
with a constant charge density n. A potential V is applied 
between both ends of the wire, and the wire is assumed 
to have a periodic structure represented by the periodic 
dielectric function e{x/a). Standard electrostatic theory 
tells us that there is a potential (j){x) such that the elec- 
tric field e{x) — —<j)'{x). The microscopic constitutive re- 
lation is d{x) = e{x/a)e{x), and the electrostatic equa- 
tion d'{x) = n gives us the equation {e{x / a)<f>' {x))' — —n, 
with the boundary conditions 0(0) — 0, (j){L) = V. For 
simplicity, we assume that L corresponds to an integer 
number of periods 27ra. This electrostatic equation has 
a unique solution (j){x) for each a (that can be obtained 
explicitly). From the explicit solution, it can be shown 
that (j){x) tends weakly to a function ^{x) when a tends 
to zero (with respect to the wire length L). The purpose 
of homogenization theory is to answer the question: to 
which differential equation is 'P{x) a solution? Again, the 
explicit solution can be used to show that <P{x) is the so- 
lution of (e^'ix))' = -n, with <P{0) = 0, <P{L) ^ V and e 
is a constant. The surprising fact is that e is not given by 
the average of e(x) over a period, but by the inverse of the 
average of l/e{x) over a period. However, such a simple re- 
lation between e(x) and e is restricted to one-dimensional 
problems, and more work is required in three dimensions. 

More generally, for a differential equation with rapidly 
oscillating coefficients, homogenization theory determines 
whether the solution has a limit when the period tends to 
zero, and to which equation the limit is a solution. This 
is exactly what we need to derive constitutive relations 
from a microscopic description of matter. Homogenization 
theory has shown that the macroscopic equation can be 
strongly different from the microscopic one. For instance, 
instantaneous micros cop ic equations can turn into equa- 
tions with memory ||2^ , ^ , local equations can develop 
non-local macroscopic terms , mixtures of optically in- 
active materials can become optically active pO|,^ 

3 Functional transformations 

In this section, some functional transformations are intro- 
duced. 

For the purpose of homogenization, a function /(r) 
is written as f(R, p), where the dependence of f is slow 
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P' 



P" 



Fig. 2. P and P' are two points in the same cell, P and P" 
are separated by a lattice vector 



for the variable R and periodic for the variable p. This 
so-called two-scale correspondence can be done explicitly 
as follows. Take a three-dimensional periodic lattice with 
Brillouin zone BZ. Let C be the Wigner-Seitz cell of the 
lattice and |C| its volume. The crystal sites will be denoted 
by Rs, and the reciprocal lattice vectors by K. 
Write the function / as a Fourier transform 

/(r) = J dqcxp(iq • r)/(q) 

= V / dqexp[z(q + K) • r]/(q + K). 

J BZ 



K 



Define now 



f(R>P)=/ rfqexp(iq • R)f(q,/9), 

J BZ 



with 



f(q,p) =^exp(zK.p)/(q + K) for q G 

K 

f (q, p) = for CiiBZ 
Then it is clear that /(r) = f(r,r), so that 



/(r) 



/ (iqexp(iq • r)f (q, r) 

JBZ 



(1) 



f(R, p) has the lattice periodicity in p and depends 
more slowly on R than on p. More precisely, homogeniza- 
tion is useful when the function /(r) varies slowly from 
P to P" (see Fig. 2), i.e. for two points that differ by a 
(short) lattice vector and varies arbitrarily from P to P' 
inside a cell. Then /(r) is transformed into f(R, p), the p 
variable describes the fast variation of the function inside 



the cell, and the R variable its smooth variation from cell 
to cell. 

This representation of a function of one variable by a 
function of a fast periodic variable p and a slower variable 
R is the essence of the two scale analysis of homogeniza- 
tion theory. 

We shall not homogenize the Maxwell equations with 
the usual two-scale function f(R, p), but with its Fourier 
transform with respect to the slow variable R: f(q, p). To 
establish a direct link between f (q, p) and /(r) and for 
future reference, we shall make use of the standard rela- 
tions for infinite Born-von Karman boundary conditions 
[ p2[ (which are derived from the Poisson summation for- 
mula II): 



\c\ 



I (iqexp(iq • R,) = 

J BZ 

I dpexp(zK • p) = |C|5k,o 
Jc 

^exp(zK.p) = |C|^% + R,) 



K 



exp(zq . R,) = ^ 5(q - K). (2) 



\C\ 



K 



An infinite Born-von Karman boundary condition avoids 
the subtle problems linked with the use of a finite Born- 
von Karman domain, for instance the question whether or 
not the other domains contribute to the reaction field. 
Using the definition of the Fourier transform we find 



f(q,p) -^exp(*K.p)/(q + K) 

K 

= (2^X1*^^P(*^'/') / c?rexp[-i(q-|-K) •r]/(r) 



K 



(2^) 
\C\ 

\C\ 

{27TY 



^ J fir/(r)exp(-iq-r)^exp[iK- (p-r)] 
X! / drf{r) exp(-iq ■ r)S{p - r + R^) 
^/(p + R,)expHq-(p + R,)]. (3) 



For a function /(r), with two-scale transform f(R, p), 
we define the average over a unit cell by 

(/(r)) = (f(R,p)) = ^y' dpf(R,p) 



\C\jc 
f^q/(q)exp(iq-R). 



BZ 



Similarly, for f (q, p) 



(f(q,p)) ^ J^dpf{ci,p) 
= /(q)- 
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Since q is in the first Brillouin zone, the average over 
a unit cell has the effect of stripping all the high fre- 
quency components off the Fourier transform of /. We 
shall use this average to go from microscopic to macro- 
scopic fields. The macroscopic fields will be microscopic 
fields averaged over a unit cell, so that all Fourier com- 
ponents of the macroscopic fields are zero when the argu- 
ment is outside the first Brillouin zone. This definition 
of macroscopic fields is common in the physical litera- 
ture 0,|ll|,|l|,|34[|3|,||l. Notice also that (f(R,p)) and 
(f (q, p)) are Fourier transforms of one another. Therefore, 
the macroscopic fields E(R) and E(q) obtained by aver- 
aging the corresponding microscopic fields remain Fourier 
transforms of one another. 

By direct substitution, it can be shown that the two- 
scale transforms of the gradient of / are (Vr -I- Vp)f (R, p) 
and («q-F Vp)f(q,p). 



4 The microscopic Maxwell equations 

To simplify the presentation, we consider a non magnetic 
sample and we neglect the effect of spin (see Ref. for the 
full theory). Moreover, the electromagnetic charges, cur- 
rents and fields have a time dependence exp(— icji) which 
will be implicit for notational convenience. 

The random phase approximation of quantum electro- 
dynamics corresponds to the following picture |Q. An 
external electromagnetic wave E°''*(r), B°''*(r) polarizes 
the dielectric crystal, creating a current density j(r) and a 
charge density 7i(r). These current and charge induce an 
electric field e(r) — E''''*(r) and a magnetic field b(r) — 
B™*(r), that induce additional current and charge, etc. 
When the medium and the field reach equilibrium, the 
vacuum (rationalized SI) Maxwell equations describe the 
connection between the total fields and the induced charge 
and current densities: 

V • e(r) = n(r)/eo 
V X e(r) = iijjh{r) 

V ■ b(r) = 

V X b(r) + itjeoMoe(r) = Aioj(i"), (4) 

where charge conservation implies V • j(r) = iuin{T). 

Linear response theory gives us the microscopic rela- 
tion 



j(r) = -iuj J dr'x{r,r') •e(r'), 



(5) 



where the electric susceptibility Xij the well-known 
expression |,|3|3|,||l 



Xy(r,r') 



-no(r)%(5(r ■ 



V;°"(r)V,"0(r') lA"0(r)t/0"(r') 



where no(r) is the charge density in the ground state, e is 
the (negative) electron charge and m its mass. The veloc- 
ity matrix elements are 

V""(r) = -^[rF*ir)VMr) - ^„(r)V.^o*(r)], (7) 

where tf'o(r) and !fVi(r) are eigenstates of the crystal. In the 
one-electron case, the sum is carried out over the occupied 
(0) and unoccupied (n) states. 

In expression (|), the first term, called the diamagnetic 
term, depends only on the electronic density, and the sec- 
ond (paramagnetic) term is usually much larger than the 
first in the optical range 

Formula ^ for the susceptibility x holds at zero tem- 
perature. The corresponding formula for finite tempera- 
ture involves a weighted sum over initial states Q. All 
the results of the present paper can be straightforwardly 
adapted to the finite temperature case. 

From the microscopic current j (r) it is customary in di- 
electric theory to define the polarization p(r) — (i/cL')j(r) 
and the displacement field d(r) = eoe(r) -I- p(r). In terms 
of those fields, the microscopic Maxwell equations become 

V X e(r) = iwb(r) 

V • b(r) = 

V • d(r) = 

V X b(r) = -icj^od(r), (8) 

where the microscopic constitutive relation is now 

d^{r) = eoe,(r) +^ / dr'xjj (r, r')ej (r'). 

j 

The microscopic fields e(r), b(r) and p(r) are, as the 
susceptibility Xij(rjr'), rapidly oscillating functions of r. 
The macroscopic fields E(R), B(R) and P(R) are aver- 
ages of the corresponding quantities over a unit cell. The 
constitutive relations problem is how to go from the re- 
lation between p(r) and e(r) to a relation between P(R) 
and E(R). Homogenization theory is a very convenient 
way to solve this problem. 



5 Homogenization of the microscopic 
Maxwell equations 

The derivation of the macroscopic Maxwel l eq uations from 
the microscopic ones is standard il^, ^ , |40| . Hence, the 
purpose of this section is just to show that homogeniza- 
tion theory gives the known results and to derive equations 
that will be used in the next sections. The homogeniza- 
tion of the Maxwell equation has been thoroughly studied 



applied mathematicians 
|, ^ and physicists pC|,|3^ 



. However, 



tl 



the microscopic 



?i(w„o - ^) 



, (6) 



constitutive relations that they used were always local: 
Xij(r,r') — Xij(r)5(r — r'). We do not need the full appa- 
ratus of homogenization theory of non-local equations [ 48| 
and we homogenize the Maxwell equations by a simple 
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adaptation of the method used for local dielectric func- 
tions. 

In this section, we follow the very clear homogcnization 
procedure of Sanchez-Palencia except that we work 
with f (q, p) instead of f (R, p). 

Using the fact that the two-scale transform of Wf{r) 
is (iq+ Vp)f(q, p), we can make the two-scale transform 
of the Maxwell equations (||) to obtain: 



iqx e(q,p) + V 
iq-b(q,p) + Vp 

iq-d(q,p) + Vp 

iq X b(q, p) 



P X e(q,p) = iujh{ci,p) 



b(q,p) = 
d(q,p) -0 



Vp X b(q,p) 



-iwpod(q,p). (9) 



For notational convenience, we do not write the fields 
e(q, p), etc, with a tilde. 

Because they vary slowly, the eventual "external" cur- 
rents and charges would appear as j(q) and n(q), with no 
dependence on p. 

Consider a scattering problem where a plane electro- 
magnetic wave is shined on the dielectric. In equations 
the order of magnitude of q will be 27r/A, where A is the 
wavelength of the incident wave. The order of magnitude 
of the unit cell dimensions is I = \C\^^^. Let a = l/X, the 
order of magnitude of Vp • e(q, p) is 1/a times the order of 
magnitude of q • e(q, p). Now we expand all fields as sums 
of the type 



e(q, p) = e(0) (q, p) + ae^' (q, p) + a^e^-''^ (q, p) + ■ ■ -(10) 

where all terms of the expansion are periodic in p. Since 
a is small, we keep only the first term of the expansion to 
define the macroscopic fields as E(q) = (e^"^ (q, p)), etc. 
The validity |^ and asymptotic convergence of this ex- 
pansion is the main technical difficulty that was solved by 
the mathematicians who homogenized the Maxwell equa- 
tions (see, for instance, Ref. [p4[ where the general term 
of expansion (|l^) is given). To be complete, we need to 
know that the order of magnitude of uj is 27rc/A, the order 
of manitudc of b is e/c and the order of magnitude of d 
is toe. 

Introducing the expansions ® into Eq.(|) and gath- 
ering all terms of same power of a we obtain four equations 
for the term: 



Vpxe(°)(q,p) =0 
Vp-b(0)(q,p) =0 
Vp.dW(q,p) =0 

Vp xbW(q,p) =0, 

and four equations for the a" term: 



(11) 



iq X e(°) (q, p) 4^ aVp x e^^) (q, p) = tuh^"^ (q, p) 
*q-b(0)(q,p)+aVp-b(i)(q,p) =0 (12) 
zq-dW(q,p)+aVp-d(i)(q,p) =0 

iq X bW(q,p) +aVp x h^'^q, p) - -zc^pod^^Hq, P), 



(recall that aVp has the same order of magnitude as q). 
These two sets of equations are sufficient to determine the 
macroscopic Maxwell equations. Further terms would be 
required if the wavelength were not very large with respect 
to the unit cell dimensions. 

To obtain the macroscopic Maxwell equations, we must 
eliminate the terms of order one, e.g. b^^)(q, p). This is 
achieved by transforming the cell average of the divergence 
into an integral over the surface dC of the unit cell: (Vp • 

b'^Hq.P)) = ^/\C\jQc'^''^^^H^yP)'^^- Oil opposite sides 
of cell C, the outgoing normal h is reversed while b^^^ (q, p) 
is equal (by periodicity). Therefore, the surface integral is 
zero and (Vp ■ b(i'(q, p)} = 0. A similar reasoning leads 
to (Vp xb(l)(q,p)) =0. 

Thus, averaging Eq. (|l^) gives = iq • (b(*'^(q, p)) + 
a(Vp-b(i)(q, p)) = jq-B(q), because the macroscopic field 
B(q) was defined as the average of h^^\q,p). Carrying 
out similar calculations for the other microscopic Maxwell 
equations we obtain the macroscopic Maxwell equations 
in momentum space: 

iq x E(q) = ia;B(q) 
^q•B(q) =0 
iq-D(q) =0 

iq X B(q) = -icjpoD(q). (13) 

Fourier transforming back to real space we obtain the 
usual macroscopic Maxwell equations: 

V X E(R) = icjB(R) 

V ■ B(R) = 

V • D(R) = 

V X B(R) = -iwpoD(R). 

Notice that this homogcnization procedure can be ap- 
plied directly to Eqs.(H), but the present formulation is 
required for consistency with the derivation of the consti- 
tutive relations. 

Homogcnization now proceeds classically. We start by 
homogenizing the magnetic field b*^"^ (q, p). From the equa- 
tion Vp X b('')(q,p) = we deduce that there is a pe- 
riodic potential 0(q, p) and a function of q, denoted by 
F(q), such that b('')(q,p) = F(q) - Vp0(q,p). Taking 
the average of both sides, and considering the periodic- 
ity of we obtain (b(°)(q,p)) = F(q) |§ and, since the 
left-hand side is defined as the macroscopic magnetic field 
B(q), we get b(°)(q,p) = B(q) - Vp0(q,p). Introduc- 
ing this equality into Vp ■ b'^'^^ (q, p) = we obtain the 
equation Z\p0(q, p) = 0. The only periodic solution of the 
Laplace equation is a constant p^ , so 4>{q, p) is a constant 
and b^°)(q, p) — B(q). In other words, the magnetic field 
does not need to be averaged, the zero-th order term of 
b(q, p) is smooth and does not depend on p. The periodic 
modulation is only reached at the next and smaller term 
ab(i)(q, p). 

Such an automatic averaging is not possible for the 
electric field. Starting from Vp x e(°^(q, p) — 0, the only 
conclusion that we can reach at this level is that there is 
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a periodic potential 0(q, p) such that e("^(q, p) = E(q) — 
Vp0(q, p). The determination of (j) will be the purpose of 
the next sections, but this result has already an interest- 
ing meaning. The zero-th term of the expansion of e(q, p) 
is not smooth. In mathematical terms, the limit of e(q, p) 
as a —> is still an oscillating function, contrary to our ex- 
ample f{x/a) = sin(a;/a) -I- h. This concept of "oscillating 
limit" was introduced b y th e Cameroonese mathematician 
G. Nguetseng in 198 9 pl| and has deeply simplified ho- 
mogenization theory ||52{. 



6 The constitutive relation 

Still, to give a complete description of the macroscopic 
electromagnetic properties of matter, we have to establish 
a correspondence between D(q) and E(q). 

To do this, we use the periodicity of for every lattice 
vector R, (r R, r' R) = Xii(i", i"') §31 ■ Notice that 
the same lattice vector must be added to both arguments 
of X- 

We start from the relation between the polarization 
and the electric field 



p(r)./dr'x(r,r').e(r'). 



Then, we apply the two-scale transformation Eq.(^ to 
find 



p(q,p) 



\c\ 



^exp[-iq- (p-l-Rs)] 



(27r)3 

j dr'x(p + R.,r')-e(r') 



Now the integral over all space is split into integrals over 
translated unit cells: 



p(q, p) 



\C\ 



^exp[-iq- (p-l-Rs)] 



(27r)3 

/ dp'x(p + R.,p' + R.') •e(p' + R,0. 
Jc 



If we introduce the two-scale expression for the electric 
field e(p) = Jg^ dqexp(iq • p)e(q, p) and use the period- 
icity of e(q, p) in p we get 

P(q> P) = exp[-iq • (p + Rs)] / dq' / dp' 

(27r)'^ ^ Jbz Jc 

exp[iq' • (p' + R,0]x(p + R., p' + R.') ■ e(q', p')- 

If we replace Rs by Rt -I- Rs' and use the periodicity of x 
we obtain 

p(q, p) = yft exp[-«q • (p + R*)] / rfq' 

(.^""j f JBZ 

/ dp' exp(iq' • p') y] exp[i(q' - q) ■ Ry] 
Jc 

x(p-t-Rt,p')-e(q',p'). 



The sum over s' is carried out with Eq.(|^) and, since q 
and q' belong to the first Brillouin zone, only the term 
K = contributes. We reach finally 

P(q, P) ^P' exp[-iq • (p - p' -t- Rs)] 

x(p-t-R.,p') -elq, p') 
= {x{p,p';(\) •e(q,p'))p', 

where the index p' designates cell average over variable p' , 
and where we have defined the two-scale transform of x 
as 

X{P: p'; q) = \C\ ^ exp[-iq • (p - p' + Rs)]x(/0 + Rs, p')- 

S 

(14) 

It can be checked that x{p^p'\^) has the lattice peri- 
odicity for each variable p and p' independently. A related 
definition was used by Ehrenreich |p4| . 

All quantities are now periodic and diagonal in q: they 
are in a suitable form for homogenization. From the rela- 
tion between polarization and electric field, we deduce the 
relation between displacement and electric field which will 
be our starting point: 

d,(q,p) = eoe,(q,p) +^(xy(p,p';q)ej(q,p'))p'. (15) 



The macroscopic constitutive relation is obtained by 
restricting all fields to the first term of expansion (^0|). If 
we do this in Eq.(|l5|) and write the electric field as the 
sum e^") (q, p) = E(q) — Vp0(q, p) we obtain 

4"^(q,P) = toE,{q)+Y{xtj{p,p'\q))p'Ej{q) - 

j 

eo9pXq,p) -^(x»j(p,p';q)9p^0(q,p'))p'- (16) 



If Eq. ( |16D is averaged over p, the fact that {dp. (j){q, p)} = 
[ pO| leads to the first step of the macroscopic constitutive 
relation: 

A(q) = eoE,{q) +Y{xtj{p, p';q))pp'Ej{q) 

3 

" Y{x^j (p, p'; q)5p^ 0(q, p'))pp' • (i7) 



Following Wiser [gTj, a connexion with the classical 
approach is possible through the introduction of a macro- 
scopic local field Eioc(q) defined by 



ixip, p'; q))pp' • Eioc(q) = (x(p, p'; q) • e(q, p')) 



pp ' 



so that D(q) = eoE(q) + (x(p, p'; q))pp' • Eioc(q). 

The next step is the determination of 0(q, p'). We take 
the microscopic equation Vp • d*^"^ (q, p) = derived in the 
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previous section (Eq.(|ll|)), and we apply it to Eq.(16). We 
obtain an equation for (/)(q, p): 

J2{^p^x^Ap^p'■^<^))p'EM)■ (18) 

This so-called cell equation determines a unique periodic 
solution (/)(q, p) with zero average over a period. 

We follow first the standard argument of homogeniza- 
tion theory. Let the three functions ^fe(q, p) {k = x, y, z) 
be the solutions of the cell equation for an electric field 
E equal to a unit vector in the direction k. Then, for 
a general electric field E(q), the potential is (/)(q, p) = 
A(q,p)-E(q). 

If we introduce this expression for 0(q, p) into Eq.(O), 
we obtain the macroscopic constitutive relation Di{q) = 
J2j ^^j{^)Ej{q) where the macroscopic dielectric constant 
€ij (q) is given by 



ey(q) = eo6ij +^{xikip,p';n)[6kj - dp'^Aj{q,p')]) 



pp' ■ 



This equation shows that the macroscopic dielectric 
constant is the average of the microscopic dielectric func- 
tion plus a correction term. In the case when the dielec- 
tric function is a sum of constant factors localized at all 
sites of a cubic lattice, it is standard exercice of homog- 
enization theory to show that the macroscopic constitu- 
tive relation becomes the Clausius-Mossotti equation (see 
Ref . [2^ , p. 45). Therefore, the correction term can be quite 
large and should not be neglected. 

Homogenization usually stops here, and the integro- 
differential cell equation dlq) can only be solved numeri- 
cally, which seems to be a reasonable task. For instance, 
the FLAPW approach could be used, where space is cut 
into non-overlapping spheres plus an interstitial region. 
</)(q, p) is then expanded over spherical harmonics (with 
suitable radial functions) in the spheres and over plane 
waves exp(zK • p) in the interstitial region. This reduces 
the cell equation to a matrix equation. 

However, the particular structure of Xij(r, r') can be 
used to give an explicit solution for (/)(q, p). This is the 
purpose of the next section. 



7 The cell problem 

The cell equation (O) can be interpreted as follows. If 
a constant electric field E(q) is applied to the dielectric, 
linear response theory tells us that it induces a periodic 
current given by Eq. (||) 

j(q,p) = -i^{x{p.p'\(\))p' ■E(q) 
and the corresponding periodic charge is 

«(q,p) = -(V^^)Vp • j(q,p) 

- -Vp- (x(p,p';q))p' •E(q). 



The last term is minus the right-hand side of Eq.(|18|). 
This periodic charge induces an additional electric field, 
that creates an additional polarization, and 4>{(\, p) is the 
periodic potential (with zero average) which represents 
the local electric field reached at equilibrium under the 
influence of the external field E(q). In other words, the cell 
equation (18) is the electrostatic equation for the potential 
created in a unit cell of the dielectric by an external field 
E(q). 

With this picture in mind, we can solve the cell equa- 
tion iteratively. We need the periodic electrostatic Green 
function G^(p) which is a solution of €qAG'^{p) = ~5{p) 
in a unit cell. Some properties of G^{p) are discussed in 
Rcf.p, p.121. 



g*{p)^y: 



K 



E 



exp(jK • p) 

eo|C|(|KP_z,y) 



47reo|p + Rs 



(19) 



where rj is an infinitesimal positive real. 

If the Green function is applied to both sides of Eq.(18) 
we obtain 

0(q, P)--^\j drdp'G* (p - T)Vr ■ x{t, p'; q) • E(q) 

+ p J dTdp'G*{p-T)S/r-x{r,p';(l)-Vp,cj,{q,p'). 

An iterative solution of this equation can be written, 
in simplified notation 

= - pG#V • X • E - j^G*V ■ X ■ VG#V • x • E + • • ■ 

If we introduce our iterative solution in Eq. ( pT[ ) we obtain 

D = eoE+(x+pX-VG#V-x + ■••)•£. (20) 

To be more explicit, we introduce the macroscopic sus- 
ceptibility Xii(q) and write Eq.(|20|) as D(q) = eoE(q) + 
x(q) • E(q), where 



x*i(q) 



CxC 



dpdp'xi]{p,p';(\) 



+ 1743 E / dpdp'dTdT'xim{p,T;q) 
dr^G*{T - r')a,;Xn,(T',p';q) + • • • (21) 



>(q) = (x(q)-[l-pVG#V.x(q)]-^)pp'. 



To sum the right-hand side of Eq. (|20|), we need a sepa- 
rable form for x, akin to work of Cho |p3|| or the coupled- 
antenna theory of Keller Q . If we neglect the diamagnetic 
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component of x in El- (j^) i we can write it in the separable 
form 



X^3 



(p,p';q) = ^/r(p,q).g;(p',q). 



(22) 



Introducing this representation in Eq. ( |20| ) , everything de- 
couples and we obtain Di{q^ = eij{q)Ej{q) where the 
macroscopic dielectric constant ey (q) is now given by 



ey(q) = eoSi 



E 

nn' 



(/r(p,q))[l-M(q)]„:,(.9/(p,q)), 



(23) 



and where the reaction field matrix (or screening matrix 
Q) is defined as 

M(q)„„, = \C\y(gnp,q)dp,G*{p-p')d,,ff{p',q)) 



PP' ■ 



(24) 



Equations ( p3| ) and (|2^), giving an explicit expression 
for the macroscopic dielectric constant, are the main result 
of the paper. These equations can be useful when only a 
few states n contribute to the susceptibility x(p, p'; q) 

Integrating Eq. (|2^ ) by parts and using the periodicity 
of all functions involved, we can derive an alternative form 
for the reaction field matrix M(q)„„/ , which decreases the 
singularity of the derivative of G"^ : 



^C\{G*{p - p')\7, ■ g"(p,q)Vp, • r'(p', q)), 



An alternative summation of the series ( B 0| ) can be 
obtained by separating G'^{p — p') with Eq.(|19D. A calcu- 
lation similar to the foregoing one leads to 



e»j(q) = £0% + {Xtj{p,p']ai))pp- 
1 



-E' 

K,K' 



eo|K||K' 



mn 



exp(zK • p')dp'^Ximip, p'; q)) 



PP 



(1 - ^)kk' (exp(-iK' • p)dp^Xnj (p, p'-. q))pp' , 
where the reaction field matrix is now 
1 



KK 



■E 



(exp(— iK • p) 



eo|K||K'| 
[9p,dp'Xij{p, p'; q)] exp(iK' • p))pp' 



(25) 
(26) 



The notation means that the sum is over all non- 
zero reciprocal lattice vectors. To show that the terms 
K = or K' = do no contribute, we reintroduce the 
infinitesimal number —irj of Eq.(|l9|) and, for instance, the 
K' = term gives us {i/v){^p'X{PiP^ ^)pp' = Oj because 
the average of a divergence is zero pQ] . Similarly, A'kk' is 
zero for K = or K' = 0. The last expression for e(q) is 
computationnaly effective when xIPtP't^) is smooth and 
only a few K contribute. 

Notice that Eq. (|2^) amounts to using the Fourier trans- 
form of the susceptibility Xij(p, p';q)- We do not discuss 
this approach further, since it has been used by many au- 
thors ||l|,|l],[l|,||,|6[. It can be checked that Eq.(||) is 
equivalent their results, in the limit where a —>■ 0. 



8 Basic properties 

In this section, some basic consequences of the macro- 
scopic constitutive relations are derived. Firstly, we show 
how the usual concept of electric dipole transition is re- 
covered, then that the macroscopic constitutive relations 
are causal and have the required symmetry properties. Fi- 
nally we prove that the average of the microscopic energy 
density is given by the macroscopic energy density and we 
discuss the reaction field matrix. 



8.1 Electric dipole transitions 

To have a constitutive relation in the real space, we back 
Fourier transform the equation D(q) = e(q) •E(q) and we 
obtain 



(27) 



where 



- !■') = -^^^ j dqexp[iq • (r - r')]e,j(q). (28) 

Equation (^7|) is typical of a homogeneous (but gen- 
erally anisotropic) medium The fact that Eq.(27) 
is non-local corresponds to spatial dispersion, which has 
been much studied by the Russian scho ol p7|l , and has 
received renewed interest recently [|57|, ^|59t|60| . Besides, 
the q-dependence of eij(q) can be observed experimen- 
tally by inelastic electron, x-ray or neutron scattering [ |6T[ 
In Eq.(p8[), it is not necessary to restrict the inte- 
gral to the first Brillouin zone since, by definition, eij(q) 
is zero outside this zone. 

If eij(q) is smooth near q = 0, it can be expanded in 
a Taylor series 

eij(q) = (0) + ^ qkdk£ij{0) + ^ qkqidkdieij{0) H 



The first term corresponds to the electric dipole ap- 
proximation, the second term describes optical activity, 
the following term corresponds to electric quadrupole and 
magnetic dipole transitions. When only the first term is 
kept, then Eq.(28) becomes (r— r') = (0)5(r— r') and 
the constitutive relation is now local D(r) = e(0) • E(r). 



8.2 Causality 

A susceptibility is causal if it satisfies the Kramers-Kronig 
relations, that we write in the following form B0| 



1 



+ 00 



27ri ./_oo uj' — Lo — ie 



-duj'. 



(29) 



Because of the structure (||), the microscopic suscepti- 
bility x(r,r';w) is causal. Thus, by hnearity, its two-scale 
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transform yfp, q' \ q, lo) is causal. To be valid, the macro- 
scopic susceptibility ( P3|) has also to be causal. 

Starting from Eq.(p9|), it can be shown by recurrence 
that, for any integer n > 0, 



1 



+ 00 



-dJ. 



LJ' — U! — te 



(30) 



If we write Eq. (Elh in simplified notation as 



Xiio) = ixicj)) + i^^ixH ■ VG#V • x{iu)) + • . • , (31) 

we can apply the Kramers-Kronig transform to the right- 
hand side of Eq.(|3l|). Then, Eq.(po[) shows that each term 
of the right-hand side is transformed into itself by the 
Kramers-Kronig transform. Therefore 



+ 00 



-duj' 



oj' — uj — le 



and the macroscopic constitutive relation is causal. 

Dolgov and coll. |64 have pointed out that x(q, cj) 
might be non-causal for a; = 0. 



D is a symmetry operation of the Bravais lattice, so that 
D^^Rs is a lattice vector and we obtain 



The final step is to replace all the x in th e rig ht-hand 
side of Eq.(^ll) by the right-hand side of Eq.(|3^. For all 
integrals we make the change of variable p = T~^p, then 
we use the fact that G*{Tp - Tp') = G*{p~ p')> the 
identity dp. /dp = Dy , the orthogonality of D and the 

—3 

periodicity of all functions involved to prove that 



(34) 



For the dielectric constant £^ (0), we restrict Eq.(34) 
to the electric dipole contribution q = and we obtain 
the expected result 



ey(0) 



jj'^i'j' 



(0) 



8.3 Symmetry 

To assert the validity of our result, it is important to check 
that the symmetry group of the macroscopic dielectric 
constant is the point group of the crystal space group. 

Because of the periodicity of x(r, r'), the lattice trans- 
lations do not intervene. The other symmetry operations 
of the crystal are T = D + 1, where _D is a (possibly im- 
proper) rotation and t is a translation shorter than the 
lattice vectors. Let T be a crystal symmetry operation, 
there is a basis of eigenstates !fVi(r) such that 

<f'„(r) =r*'„(r) =if„(r-ir). 

Using this property we obtain the following transforma- 
tion rule for x(r, r') 

X^Ar,r') ^Y.^^^'^n'Xer{T-\T-h'). (32) 
i'j' 

The presence of the matrices Dai comes from the deriva- 
tives in definition (|^) , which get rid of the translation part 
of T. 

Using Eq.dl) in Eq.(|l|), we find 
X«j(p,p';q) = |C|^exp[-iq- (p-p' + R^)] 



We make the transformation T~^{p + Rg) = D~^{p + 
Us) — t = T~^p -\- £)~^Rs, we write the argument of the 
exponential function as q • (p — p' -I- R^) — Dc\ ■ {T^^p — 
T~^p' +D~^'Eis), we use the fact that, for any space group. 



8.4 Energy density 

The microscopic energy balance is ^t\ : 
- J drj(r, t) ■ e(r, ^ J drh{r, t) 



dh{v,t) 
dt 



-eo f dre(r, t) ■ ^^^^ + — / ^(i-, t) x b(r, t) ■ da. 
J at po 



The definition of d(r) gives us 

5e(r,i) ad(r,t) 



eo- 



j(r,t), 



dt dt 
and the energy balance can be written 

j dru{r, t)+ J s(r, t) ■ da = 0, 

where the microscopic electromagnetic energy is 

and the microscopic Poynting vector is 

s{r,t) — —e{r.t) x b(r,t). 
Po 

If we restrict the definition of the macroscopic fields to 
the average of the first term in expansion ( p^ ) , we want to 
determine whether the average of the microscopic Poynt- 
ing vector is equal to the macroscopic Poynting vector 
S(r, i) = ^E(r, i) x B(r, i), and whether the average of 
the microscopic energy density is equal to the macroscopic 
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energy density. If this were not the case, the energy ar- 
guments using the macroscopic Maxwell equations would 
lack any microscopic basis. 

For the Poynting vector, the answer is immediately yes, 
because we have shown that b^''^(q, p) = B(q). Therefore 
b(°)(r, = B(r,f), the magnetic field does not oscillate 
rapidly, it is equal to its average and 

(s(r,t)) = -(e(")(r,<)xb(°)(r,i)) 
Mo 

^ J_(e(")(r,<)} X B(r,t) = — E(r,i) x B(r,t) 
Mo Mo 

-S(r,t). 

The proof is similar for the magnetic part of the en- 
ergy density. However, the problem is more difficult for 
the electric part, because it is a product of two functions 
that oscillate very rapidly, and it is not obvious that the 
average of the product is equal to the product of the av- 
erages. We prove now that the average of the microscopic 
energy density is indeed given by the usual macroscopic 
formula. 

First, we Fourier transform the fields in space and time 

as 



8.5 Reaction field matrix 

It is also possible to integrate by parts to apply both 
gradients to the Green function, but the double gradi- 
ent of the Green function has a singularity that must be 
treated with care [^|7^. On the other hand, this form 
has the advantage of recovering the usual dipole-dipole 
interaction of classical dielectric theory as will be 
shown now. In Eq.(|2j), an integration by parts trans- 
fers the derivative 9p' from /" (p', q) to G*(m — p'), then 
dpi G'^{p — p') = —dp.G'^{p — p'). Therefore, if we define 
the matrix Green function G'^{p — p') by 

G*{p-p')=dp^dp^G*{p-p'), 

The definition of the reaction field matrix becomes 

A./(q)„„, = |C|5](g"(p,q)-G#(p-p')-f"'(M',q))pp'- 



To proceed, we write 

G*{p p') = G"(p - m') + E -P' + ^S 



.(0) 



(r,i) = 



/ dq [ duj exp[i(q • r - ujt)]e^"'> (q, r) Weiglhofer has shown that [|9[|7C 

Jbz J 



(the variable uj is still implicit for the fields). 

We need to show that (e(°)(q',p) • d('')(q, p)} is equal 
to the corresponding macroscopic product E(q') • D(q). 

The microscopic displacement field d'^°^(r, q) is given 
by Eq.(|l|) and the eletric field by e(o)(q,p) = E(q) - 
Vp0(ci, p). The macroscopic displacement field is given by 
Eq.(|l3), so that 

(e(0)(q',p).d(0)(q,p)) =E(q')-D(q) 

-(Vp(^(q',p)-x(M,M';q))pp' -ECq) 

+eo(Vp0(q',p)- Vp0(q, p)} 

+ (Vp</>(q',p) • xip,p';q) ■ Vp'(/)(q,p'))pp', 

where we have used the fact that (Vp0(q, p)} = 0. Inte- 
grating by parts to eliminate the gradient of (/i(q',p) we 
obtain 

(e(0)(q',p).d(")(q,p)) =E(q')-D(q) 
+ ((^(q',p)[Vp-x(M, p';q)-E(q) 
-eo^p0(q, m) - Vp • x{p, p'; q) • Vp/(/>(q, p')])pp' 
= E(q') . D(q), 

where the last step was derived using the cell equation 
(II). 

A mathematical study of the convergence of the en- 
ergy density was carried out for a local susceptibility by 
Markowich and Poupaud p5[ . 

When magnetic properties are taken into account, that 
simple result could fail. If this were the case, we would 
have a microscopic basis for the non-standard Poynting 
vectors investigated in Ref.p§]. 



G°,(r) 



-7—dr,dr.- 

47reo r 



47ren 



The first term gives the depolarization term of the classical 
Lorentz theory j4^, the second one is more delicate and 
is not considered in the classical approach. 

Finally, we shall need the multiple-scattering expres- 
sion for G°(p — p' + Rs), vahd when \p\ + \p'\ < \ILs\ (see 
Rcfs.pi^,ra for this and the most general cases). 



Go(p-p' + R.)= E Yrip)Mp)He,ni' m'iRs) 

Y/r'*ip')n,{p') 

where je{p) = p7(2^ + 1)!! and 

He^i>jn'(Rs) = -47r(-l)^G^^+^,„,_„Yg™^7'"(Rs) 



{2£ + 2£' - 1)!! 



e+i'+i 



(35) 

is a Gaunt coef- 



with (-1)!! = 1 and where G^„^+^,„,_„ 
ficient. 

Because of the derivatives in the definition of G'^{p — 
p'), the terms £ = oi £'=0 in Eq.(|35|) do not contribute. 
Therefore, the cell-to-cell electrostatic interaction begins 
with a 1/Rg. In other words, there is no interaction be- 
tween the cell total electron charges (that would be in 
1/Rs), although the cell are not neutral. To be more ex- 
plicit, the cells are neutral if one consideres the nuclei. 



Christian Brouder, Stephanie Rossano: Microscopic Calculation of the Constitutive Relations 



11 



but the nuclei charges did not enter our description, and 
their contribution to the dielectric function is very small 
at optical energies. Therefore, there is no 1/Rs interaction 
although only electrons are taken into account. 

For an octahedral or a tetrahedral symmetry, group 
theory tells us that Y^itls) is zero for ^=1,2 and 3. 
Therefore, the electrostatic interaction starts with a 
dependence. 



9 Bloch wave functions 

Many works are devoted to the calculation of the optical 
response of semiconductors, using a band-structure ap- 
proach. To deal with this case, we specify now our results 
to the case of one-electron wavefunctions. 

The Bloch theorem tells us that one-electron wavefunc- 
tions in a crystal can be written exp(«k • r)M;(r; k), where 
the index I runs over the bands. When this equation is 
introduced into the velocity matrix elements (m), we find 



ifi 



V" (r;k,k') = -— exp[z(k'-k).r] 
2 m 

(r^iMc' -n;,k)W"'(r;k,k'), 

where nj^k is the occupation of the l-ih band with Bloch 
vector k and the periodic functions (r; k, k') are 

W"' (r; k, k') - u*i (r; k) Vu,, (r; k') - m, (r; k') Vt/^ (r; k) 
+i(k + k')ur(r;k)ur(r;k'). 

Neglecting the first and the last term in Eq.(|^) we find 
X., (r,r') = I dkdk'exp[z(k' - k) • (r - r')] 



W^" (r;k,k')W^j '(r;k',k) 
Ei,{k')-Ei{k)-hLj 



We treat now the electric dipole case q = 0. From the 
definition of the two-scale transform of x and from identity 
(0) we find 



J 71/ 



ni'M - it-im) (36) 



wfip-Xk)wpip'-xk 



(37) 



Therefore, in the electric dipole case, the transitions 
are vertical. When spatial dispersion is investigated, then 
the transitions are no longer vertical, they explore a part 
of the bands around the vertical . 

We can use the Schrodinger equation for ui and u/' to 
show that 



In Ref.|5j, p. 121, Keller has proved that such a relation 
between V • V'^"(r) and matrix elements of the density 
operator holds also valid in the many-body case. 

To use the macroscopic constitutive relation (E3|) we 
must define a separable form ( p^ ) for x in Eq.(|37[). We 
can choose 



f- (p;k) 



(2vr) 



HI' 



-(^^;^k - ";,k) 



Wt (p;k,k) 



Ei,(k)-Ei{k)-nij' 
fff (p';k) = W^'V;k,k). 

Finally, using identity (|3^) and the alternative definition 
of the reaction field matrix M we find 



kk 



\C\{2n 



{ni' 



(E, jk) - Eiik)){Ei_{k) - Eij jk)) 
{Ei,{k)~Ei{k)~hu;) 
{u*i{p;k)ui>{p;k)G*{p - p)u}{p' ■,'\^ui_{p ■,'^) . 

Because of the presence of the Bloch vector indices k and 
k, the reaction field matrix is huge and its inversion will 
be a heavy computation. 

All the ingredients are now given for a band-structure 
calculation of the macroscopic dielectric constant in the 
electric dipole approximation. The final formula is 



dkdk [nil 

ll'lV 



"-i,k) 



(W,"'(p;k,k))[(l 



M)-tk--Vr(p;k,k)) 



{Ei,{k)-Ei{k)-nuo) 

A similar formula can be obtained for spatial disper- 
sion (q 7^ 0), adapting the implementation described in 
Ref . 173] , where the Coulomb singularity and the Umklapp 
processes are treated in detail. 



10 Many-body dielectric function 

Before calculating the dielectric constant corresponding to 
a many-body susceptibility, we must show that the trans- 
verse part of the microscopic electric field varies slowly. 
All fields can be written as the sum of a transverse and 
a longitudinal components. The transverse component of 
e(q, r) is defined by V-eT(r) = and VxeT(r) = Vxe(r). 
The corresponding equations for the two-scale transforms 
are 

iq • eT(q, p) -f Vp • eT(q, p) = 

(iq + Vp) X eT(q, p) = (iq + Vp) x e(q, p). 

Using the expansion ([l^) for the full and transverse 
electric fields, we obtain for the a^^ term; 



Vp • W"' {p- k, k) - - |J [El, (k) - El {k))u* {p- k)ui, (p; k) . 
Ti 

(38) 



Vp-e^°)(q,p) -0 
Vpxe^°)(q,p)-Vpxe(o)(q,p)=0 



(39) 
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where the last equality was obtained from the first of 
Eqs.(|rT|). Equations ( p9[ ) are the same as the equations 
for h^^^ (q, p) in section 4. Therefore, the conclusion is the 

same and e^^q, p) = F(q) does not depend on p, it is a 
slowly varying function. 

To determine F(q) we use the following equation given 
by the a*' term 

zq • e^^ (q, p) + aV p ■ e:p (q, p) = 0. 

The average of both sides yields q • F(q) = and F(q) is 
a macroscopically transverse vector. The other equation 
given by the aP term is 

iq X e^^ (q, p) + Vp x e^^^ (q, p) = iq x e^"' (q, p) 

+Vp X e(i)(q,p)- 

Taking again the average of both sides we obtain q x 
F(q) = q X E(q), so that F(q) = ET(q), and the micro- 
scopic transverse electric field is equal to the macroscopic 
transverse electric field. 

Notice that this conclusion follows from the Maxwell 
equations and does not depend on the constitutive rela- 
tions. 

In many-body theories, [|[|7^,^,^,|77| the induced 
current is related to the external electric field or to the 
transverse part of the total electric field (see Ref.||] for 
a detailed derivation). Since the experimental dielectric 
function relates the macroscopic electric displacement to 
the macroscopic total electric field, we choose the latter 
formulation and write 

d(r) — eoe(r) + j dr'a{r,r') ■ eT{r'). 

Repeating all the steps of section 6, we arrive at 

d(q, p) = eoe(q, p) + {a{p, p'; q) • eT(q, p'))p' ■ 

Since eT(q, p') — ET(q) does not depend on p', the aver- 
age yields simply 

D(q) = eoE(q) -|- (ci(p, p'; q))pp' • ET(q). 

ET(q) can be expressed in terms of E(q) by 

E,(q) = E(q) - i^Mq 

r 

and the final formula for the many-body case is D(q) 
e(q) • E(q) with 



(q) = eoSr 



k 



,aik{p,p;q))pp' — -■ 



This formula looks simpler than for the case of the 
random phase approximation, but the calculation of the 
many-body susceptibility a(r, r') is much more difficult, 
since it must account for all electron-electron interactions, 
including the reaction fields. Notice also that the limit 
q — *■ is ambiguous because constant electric field cannot 
be uniquely decomposed into a transverse and a longitu- 
dinal parts. 



11 Conclusion 

Homogenization theory is usually used to calculate the 
properties of "real" materials (porous, fibrous, disordered, 
etc.). It is a somewhat just reward that applied physics can 
also be useful to basic physics. 

Here, homogenization theory was used to calculate the 
macroscopic dielectric constant from the microscopic di- 
electric function. Compared to previous works, our ap- 
proach is not restricted to cubic materials and provides 
new equations to describe the local field effect in dielectrics. 
This article does not exhaust the prospects of homogeniza- 
tion theory in solid-state physics. It is rather a detailed 
presentation of the simplest possible case and the present 
study can be developed in many directions. 

A complete derivation of the macroscopic Maxwell equa- 
tions should also take magnetism into account. In par- 
ticular, the mysterious relation between microscopic and 
macroscopic magnetic properties [7^ could be handled 
with homogenization theory, as well as the question of the 
general form of the constitutive relations in bianisotropic 
media p9| . 

Further terms of the expansion can be calculated [|4| 
to investigate the case when the wavelength of the incident 
wave is not very large as compared to the unit cell, as in 
the VUV range and for some near-field optics or inelastic 
scattering experiments. 

We have considered an infinite crystal, but homoge- 
nization theory can also treat finite crystals |47||5l] ,|5^,[79|. 
This is particularly interesting when the medium exhibits 
spatial dispersion, i.e. g ^ in Eq.([2^), because addi- 
tional boundary conditions can be required to determine 
the waves inside the dielectric body |g7,53,5^. Homog- 
enization theory is well suited to describe the boundary 
layer that forms at the surface of the dielectric body, and 
to derive the corresponding boundary conditions. To do 
this, one adds to the bulk (periodic) functions a boundary 
function that decreases exponentially out of the dielectric 
[0|8^. Besides, the Bloch decomposition of the electro- 
magnetic field seems to be a promissing alternative for 
that purpose [||,|7|,|l|. 

The present work was carried out within the linear 
response approximation. Homogenization theory is fully 
developed to deal with non-hnear equations 22, 47|,p^. 

We have considered periodic media, but homogeniza- 
tion theory applies also when the periodic structure varies 
slowly or in disordered or polycrystalline media . 

Finally, the method is not restricted to electrodynam- 
ics and can be used to calculate any constitutive relation 
corresponding to a microscopic non-local equation. 
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